Minimal Bending Energies of Bilayer Polyhedra 
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Motivated by recent experiments on bilayer polyliedra composed of amphipliilic molecules, we 
study the elastic bending energies of bilayer vesicles forming polyhedral shapes. Allowing for segre- 
gation of excess amphiphiles along the ridges of polyhedra, we find that bilayer polyhedra can indeed 
have lower bending energies than spherical bilayer vesicles. However, our analysis also implies that, 
contrary to what has been suggested on the basis of experiments, the snub dodecahedron, rather 
than the icosahedron, generally represents the energetically favorable shape of bilayer polyhedra. 
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In an aqueous environment, amphiphilic molecules 
such as lipids are observed to self-organize into bilayer 
vesicles P, [Ij, thus forming the physical basis for cell 
membranes. Bilayer vesicles generally exhibit shapes 
with constant or smoothly varying curvature [3, Q ■ But 
in recent experiments [HI, Q , bilayer vesicles with poly- 
hedral shape, consisting of flat faces connected by ridges 
and vertices with high local curvature, have been ob- 
served. In these experiments, two types of oppositely 
charged, single-tailed amphiphiles were used, with a 
slight excess of one amphiphile species over the other. At 
high temperatures, the amphiphiles were found to form 
spherical bilayer vesicles. However, provided that the 
number of excess, unpaired amphiphiles was tuned to 
some optimal range, cooling the system below the chain 
melting temperature yielded polyhedral bilayer vesicles. 
It was reported that the bilayer polyhedra were stable 
over weeks, and that their shape was consistently repro- 
duced upon thermal cycling. Furthermore, it was sug- 
gested p, Q that the observed polyhedra had icosahe- 
dral symmetry, although some uncertainty regarding the 
polyhedral symmetry remained. 

What is the mechanism governing the formation and 
symmetry of bilayer polyhedra? It was argued [H, 0| 
on the basis of the experimental phenomenology that 
elastic contributions to the polyhedron free energy dom- 
inate over entropic or electrostatic contributions, and 
that minimization of elastic bending energy alone de- 
termines the shape of bilayer polyhedra. In this Let- 
ter, we take these intriguing observations as our starting 
point and address the general problem of finding poly- 
hedral shapes with minimal bending energy. Questions 
regarding the minimal energy shape of bilayer vesicles are 
commonly answered using a variational approach [ll-Q. 
However, a given polyhedral shape is defined by the geo- 
metric parameters characterizing its vertices and ridges. 
Thus, polyhedra are inherently of singular nature, which 
severely restricts the applicability of variational calculus. 
We therefore employ a complementary method, in which 
we allow for ridges and vertices with arbitrary geometric 
properties and, on this basis, calculate polyhedron bend- 
ing energies as a function of polyhedron symmetry. In 



the remainder of this Letter, we first consider the most 
straightforward case of bilayers with uniform composi- 
tion, and then turn to the richer case in which there is 
segregation of excess amphiphiles. 

The solution of the two-dimensional equations of elas- 
ticity is a formidable challenge, and has only been 
achieved for the vertices and ridges of polyhedra in cer- 
tain limiting cases [7-ll|. Thus, in order to determine 
the clastic energies of arbitrary polyhedral shapes, we 
mainly employ simple expressions based on the Helfrich- 
Canham-Evans free energy of bending [l|- 
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where K}, is the bilayer bending rigidity, Ri and R2 are 
the two principal radii of curvature, and Hq is the spon- 
taneous curvature. The resulting expressions for poly- 
hedron energies are intuitive and only involve a few pa- 
rameters but, ultimately, are purely phenomcnological. 
We assess their validity by making comparisons to poly- 
hedron energies obtained for the aforementioned limit- 
ing cases of the equations of elasticity, which allow for 
stretching as well as bending deformations. 

Figure (Ua) shows schematic illustrations of a bilayer 
bending gradually (left panel) and sharply (right panel) 
along a ridge with dihedral angle a^. The first model is 
inspired by the electron micrographs of bilayer polyhedra 
in Refs. 0,0, [Hj, while the second model provides a more 
faithful representation of the polyhedral geometry. Based 
on the picture presented in the left panel of Fig. [IJa) , we 
approximate ridges by a bilayer bending partially around 
a cylinder of radius Ri = d/{Tr — ai), where d is the 
arc length. Following the right panel, we discretize the 
bilayer with a lattice spacing b, and assume a harmonic 
potential for the angle between adjacent bond vectors. 
Upon setting d = b one finds from either approach a 
ridge energy similar to the expression used in Ref. @ , 
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where Kb ~ Kb/b, I is the ridge length, and we have 
assumed that Hq = 0. A simple expression for the bend- 
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FIG. 1: (color online). Illustration of the contributions to 
the elastic bending energies of polyhedra: (a) Side view of a 
ridge with dihedral angle ai, (b) vertex with face angle 
(c) cross section of half of a pore around the tip of a cone (see 
inset) with apex angle n — 29 and radius r, (d) top-down view 
(left panel) and side view (right panel) of a pore composed of 
straight edges along each face. 

ing energy associated with closed bilayer vertices [see 
Fig. [TJb)] is obtained from Eq. ^ following analogous 
steps, leading to the vertex energy 

G. = f^E(^-/3.)'. (3) 

i 

where /3j denotes the face angle subtended by two ridges 
meeting at a given polyhedron vertex. 

Closed bilayer vertices may break up to form pores 
@, @: [13 1 which was suggested jH as a mechanism for 
avoiding the curvature singularity associated with closed 
vertices. Figure [Tjcjd) show two models for pores at 
the vertices of polyhedra which, similarly as before, are 
inspired by the experimental images in Refs. [i, 0, [13 
[Fig. (He)] and a stricter interpretation of the polyhe- 
dral geometry of bilayer vesicles [Fig. [Ijd)]. In our first 
model [see Fig.[ljc)] we approximate the vertex of a given 
polyhedron by a cone with apex angle tt — 29, where 
6 = tt/2 — arccos(l — J7/27r) for a solid angle $7 sub- 
tended by the polyhedron vertex. In our second model 
we assume that, along each face, the pore consists of a 
straight cylindrical edge [see Fig.[IJd), left panel], which 
bends through an angle 7j across a ridge from one face 



to a neighboring face [see Fig. [Ijd), right panel]. In 
both cases, the elastic pore energy Gp can be evaluated 
on the basis of Eq. ([T]) and is found [3] to depend on 
the pore radius, the monolayer bending rigidity, K^, and 
the monolayer spontaneous curvature, Hq, as well as the 
face angle and solid angle characterizing the geometry of 
a given polyhedron. For physically relevant parameter 
ranges, Gp increases with decreasing Hq and fi. 

We now turn to the presence of excess amphiphiles in 
bilayer polyhedra. The expulsion of excess amphiphiles 
from flat bilayers and the resulting molecular segregation, 
together with the high spontaneous curvature of single- 
tailed excess amphiphiles, are thought (1, 0] to have two 
principal effects on the bending energies of polyhedra. 
On the one hand, excess amphiphiles can seed pores into 
bilayers 0,0] and, thus, pores may have a role beyond re- 
ducing the elastic energy of polyhedron vertices. In par- 
ticular, it has been suggested [6[ that excess amphiphiles 
produce pores in the spherical bilayer vesicles from which 
bilayer polyhedra originate upon cooling. On the other 
hand, it has been found Q that excess amphiphiles pref- 
erentially accumulate along the ridges of polyhedra. As 
a result, molecular segregation can decrease the bending 
energy of the outer monolayer at ridges. 

The above observations suggest a simple description of 
how vertex and ridge energies are modified by the pres- 
ence of excess amphiphiles. Ideally, excess amphiphiles 
are arranged along ridges such that they induce an 
anisotropic spontaneous curvature commensurate with 
the dihedral angle. Assuming such "perfect segregation" , 
only the bending energy of the inner layer must be con- 
sidered when computing ridge energies. We therefore ob- 
tain a lower bound on the modified ridge energy which 
takes a similar form as Eq. ([2]), but with the rescaled 
bilayer bending modulus Kb replaced by the rescaled 
monolayer bending modulus = K^/b. Thus, pro- 
vided that the optimal amount of excess amphiphiles is 
present [1, 0], the ridge energy is lowered by a factor 
K^/Kb, with the number of pores seeded into spherical 
vesicles equal to or greater than the number of polyhe- 
dron vertices. Experiments 0, 0] and simulations [l^ 
suggest that Kl/K^ ^ 10^^^ 

From the simple model for perfectly segregated bilayer 
polyhedra described above one finds tliat the optimal ra- 
tio of the amphiphile species in excess to the total am- 
phiphile content is given by r/ « 0.51 for the polyhedron 
sizes observed in experiments 0, 0] . This optimal value 
for r/ is a direct result of amphiphile and polyhedron 
geometry and, hence, does not depend on any elastic pa- 
rameters. Tlie corresponding experimental estimate is 
ri « 0.57 0,0]. We expect that in experiments not all 
excess amphiphiles are segregated along the ridges and 
vertices of polyhedra as a result of, for instance, entropic 
mixing within bilayer polyhedra or the formation of mi- 
celles [6|. Thus, our theoretical estimate for r/ is in broad 
agreement with experimental observations. 
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For certain limits of the equations of elasticity, approx- 
imate solutions corresponding to polyhedron vertex 0- 
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and ridge [10|, [11 



energies have been obtained. In 
particular, it has been found that 5-fold disclinations 
in hexagonal lattices are accommodated for small lat- 
tice sizes through a stretching of lattice vectors. How- 
ever, for large enough lattice sizes it becomes energet- 
ically favorable to buckle out of the plane 01 , in which 
case the energetics of the system are dominated by bend- 
ing. This behavior is characterized by a dimension- 
less quantity known as the Fbppl-von Karman number 
r = YR^ /Kb, where Y is the two-dimensional Young's 
modulus and R is the lattice size. As F — > cxj, pro- 
nounced ridges develop [1, 0] between the 5-fold disclina- 
tion sites of icosadeltahedral triangulations of the sphere. 
In this limit, the total elastic ener gy is dominated by 
ridges and was determined in Refs. jiol . [ll| to be of the 
form ifh(7r - af/^l^/^f{Y/Kb), with 0.1 <f(Y/Kb) ^ 
0.4 nm-i/3 foj. bilayer polyhedra 0, [STIli- Thus^ 
the asymptotic expression for the ridge energy found in 
Refs. [iS ini leads to a similar dependence on the dihe- 
dral angle and proportionality factor as in Eq. while 
increasing sub- linearly with the ridge length I. 

The lowest energy states of icosadeltahedral triangu- 
lations of the sphere are found to resemble icosahedra 
for F ^ 10^ 0, [^, which corresponds to a vertex en- 
ergy greater than Si^Tf, with, for instance, a value \2Kb 
for F = 10^". This compares quite favorably with the 
estimate G„ « llA'f, implied by Eq. ([3]) for the icosa- 
hedron. Moreover, we find [l^l that our two models for 
pores predict similar ranges for the pore energy Gp, with 
the competition between pores and closed bilayer vertices 
governed by the ratio K^/Kb- In particular, the afore- 
mentioned estimate K^/Kb £ lO'^ [i[i,[l3 

implies that 

closed bilayer vertices will be unstable to the formation of 
(closed) pores, thus removing the singularity associated 
with polyhedron vertices. This is consistent with exper- 
imental observations [^, H, [l^l and allows adjustment of 
the volume of bilayer polyhedra for a fixed total area or 
number of amphiphiles. 

Ridges impose an energetic cost and, hence, one ex- 
pects that for a fixed area and dihedral angle the faces of 
polyhedra relax to form regular polygons. We therefore 
focus here on the convex polyhedra with regular polygons 
as faces, but we have also considered polyhedra with ir- 
regular faces [l^. The class of convex polyhedra with 
regular faces encompasses the five Platonic solids, the 
thirteen Archimedean solids, the two (infinitely large) 
families of prisms and antiprisms, and the 92 Johnson 



solids [la ll6| . It has been shown [17| that this list ex- 
hausts all convex polyhedra with regular faces. Thus, 
counting prisms and antiprisms as one solid each, there 
are exactly 112 convex polyhedra with regular faces. 
With each of these polyhedra a specific set of parame- 
ters characterizing ridges, vertices, and faces is associ- 
ated, leading to distinct contributions to the total clastic 
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FIG. 2: (color online). Total elastic bending energies of the 
convex polyhedra with regular faces, obtained for the case of 
perfect segregation of excess amphiphiles, and bending energy 
of the sphere, normalized by the total bending energy of the 
icosahedron, d, with (a) pores with r — nm and (b) pores 
with r = 20 nm at each polyhedron vertex. The snub dodec- 
ahedron corresponds to the bold curve minimizing bending 
energy in (a,b) for Rp ~ 500 nm. 



bending energy. 

We have evaluated the total clastic bending energies 
of all convex polyhedra with regular faces. As illus- 
trated in Fig. [5J polyhedron energies are compared by 
plotting elastic energy as a function of the polyhedron 
radius Rp which is related to the polyhedron area 
A via A = 47ri?p, such that, for each value of Rp, all 
shapes have the same total area. While the quantita- 
tive details of the resulting energy curves depend on the 
particular combination of the aforementioned expressions 
for ridge, vertex, and pore energies used, we find that all 
curves share the same basic qualitative features. Con- 
sistent with a previous study 0], the icosahedron [see 
Fig. El^a)] minimizes bending energy among the Platonic 
solids. However, we also find that, in general, the icosahe- 
dron docs not minimize bending energy among arbitrary 
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FIG. 3: (color online). Image representations of (a) the icosa- 
hedron, (b) the snub dodecahedron, (c) the snub cube, and 
(d) the great rhombicosidodecahdron. The polyhedra in (b) 
and (c) are chiral. 



polyhedral shapes. In fact, for large enough polyhedron 
sizes, the snub dodecahedron [see Fig. El^b)] is the poly- 
hedral shape minimizing bending energy among the eon- 
vex polyhedra with regular faces, and the snub cube [see 
Fig. El^c)] also has a lower energy than the icosahedron 
in this limit. 

Allowing for an optimal number of excess amphiphiles, 
we find (see Fig. [5]) that polyhedra can have lower bend- 
ing energies than the sphere, but only if we permit molec- 
ular segregation along ridges as observed in Ref. Q. 
Segregation at pores, which was originally suggested in 
Ref. [5| as a potential mechanism stabilizing polyhedral 
shapes, is not sufficient to produce polyhedra with bend- 
ing energies which are favorable compared to the sphere 
[l3|. Indeed, if we assume that pores are closed, bi- 
layer polyhedra are energetically favorable for the exper- 
imentally observed polyhedron radius Rp ~ 500 nm [see 
Fig. [2Ia)], with the snub dodecahedron as the minimum 
energy shape among the convex polyhedra with regular 
faces. If we allow pores of a finite size, a sequence of poly- 
hedral shapes is obtained which minimize bending energy 
for smaller polyhedron radii [see Fig.[5Jb)]. The most no- 
table of these polyhedral shapes is the great rhombicosi- 
dodecahedron [see Fig. El^d)], which surpasses the snub 
dodecahedron in bending energy at Rp w 300 nm. How- 
ever, according to our analysis, the snub dodecahedron 
represents the minimum in bending energy among the 
convex polyhedra with regular faces for the polyhedron 
sizes and pore sizes (r w 20 nm) found experimentally 

In summary, we have used Eq. ([T]) and a variety of 
other expressions fMu\ to systematically evaluate the 
bending energies of bilayer polyhedra. We find that, con- 
trary to what has been suggested on the basis of exper- 
iments [HI, Q , the snub dodecahedron and the snub cube 



generally have lower total elastic bending energies than 
the icosahedron. This result is consistent with several 
complementary theoretical studies which suggest that the 
elastic energies of chiral shapes such as the snub dodec- 
ahedron and the snub cube can be favorable compared 
to the icosahedron [l8l-l20l| and that, even if the icosa- 
hedral shape is imposed, the minimum energy structure 
may still be chiral [2l|. While we followed here the ex- 
perimental phenomenology 0, 0] and assumed that min- 
imization of bending energy governs the shape of bilayer 
polyhedra, other contributions to the free energy, as well 
as kinetic effects [13, Ull , could, in principle, modify the 
preferred polyhedral symmetry. In light of our results, 
we suggest revisiting the symmetry of bilayer polyhedra, 
and the thermodynamic or kinetic mechanisms poten- 
tially governing their formation and stability, in greater 
experimental detail. 
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